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The present paper is devoted to the relativistic statistical theory, introduced in Phys. Rev. E 
66 (2002) 056125 and Phys. Rev. E 72 (2005) 036108, predicting the particle distribution lunction 
p{E) — exp^(— — fi]) with exp^(a::) = (\/l + k^x^ + nxy^"", and < 1. This, experimentally 
observed, relativistic distribution, at low energies behaves as the exponential, Maxwell-Boltzmann 
classical distribution, while at high energies presents power law tails. Here, we obtain the evolution 
equation, conducting asymptotically to the above distribution, by using a new deductive procedure, 
starting from the relativistic BBGKY hierarchy and by employing the relativistic molecular chaos 
hypothesis. 



I. INTRODUCTION 

In the last few decades, the power-tailed statistical dis- 
tributions have been observed in a variety of physical, 
natural or artificial systems. In [l| one can find an up- 
dated discussion of power-law tailed distributions and an 
extensive list of systems where that distributions have 
been empirically observed. 

For instance the cosmic ray spectrum, oc xiP^i ~ 
/3/i), obeys the Boltzmann law of classical statistical me- 
chanics i.e. x{^) ^ exp(— x) at low energies, while at 

a:— f 

high energies this spectrum presents power law fat tails 
i.e x(a;) ~ The cosmic rays problem was ap- 

proached for the first time, by using a different distribu- 
tion from the Boltzmann one, in 1968. In his proposal 
Vasyliunas heuristically identified the function x{^) with 
the Student distribution function which presents power- 
law tails [l|. 

Currently, there is an intense debate regarding the 
physical origin of the experimentally observed non- 
Boltzmannian distributions. Recently, after noting that 
the power-law tails are placed in the high energy region, 
where the particles are relativistic, the question has been 
posed whether the solution of the problem, i.e. the theo- 
retic determination of the function x(a;) and consequently 
of the related distribution and entropy, can be explained 
by invoking the basic principles of special relativity. Few 
years ago in refs. 0, 0] has been obtained within the 
special relativity the following simple expression for the 
function x{^)^ i-^- x{^) = ^xp^(a;) with 



exp, 



(1.1) 



In the last few years various authors have considered 
the foundations of this statistical theory, e.g. the H- 
theorem and the molecular chaos hypothesis 0, [1] , the 
thermodynamic stability 0, , the Lesche stability [iTI - 
[T^ , the Legendre structure of the ensued thermodynam- 
ics [l^ 111] , the geometrical structure of the theory [13] , 
etc. On the other hand specific applications to physi- 
cal systems have been considered, e.g . the cosmic rays 
0], relativistic [isj and classical [l9'| plasmas in pres- 
ence of external electromagnetic fields, the relaxation 
in relativistic plasmas under wave-particle interactions 
pol [2H , anomalous difFusioii_ 2^, [231 > kinetics of inter- 
acting atoms and photons [24| . particle kinetics in the 
presence of temperature gradients l25ll. particle systems 
in external conservative force fields , stellar distribu- 
tions in astrophysics [13, [l^, quark-gluon plasma for- 
mation [29| . quantum hadrodynamics models [3d] , the 
fracture propagation [sij, etc. Other applications regard 
dynamical systems at the edge of chaos [32-341, fractal 
systems [35], the random matrix theory |36ll37j. the error 
theory |38|] , the game theory ^] , the information theory 
[ioj . etc. Also applications to economic systems have 
been considered e.g. to study the personal income distri- 
bution [4ll - |43| . to model deterministic heterogeneity in 
tastes and product differentiation [3 [HI etc. 

In the present contribution we reconsider critically the 
kinetic foundations of the statistical theory based on the 
function exp^(a:). Our main goal is to show that the stan- 
dard principles of ordinary relativistic statistical physics, 
i.e. kinetic equation, H-theorem, Molecular Chaos Hy- 
pothesis (MCH) , conduct unambiguously to the relativis- 
tic generalization of the classical Boltzmann entropy and 
Maxwell-Boltzmann distribution. 



The parameter < 1 is the reciprocal of light speed, in 
a dimensionless form, while 1/k^ is the particle rest en- 
ergy. In the classical limit k — > 0, exp^(a;) reduces to the 
ordinary exponential function, defining the Boltzmann 
factor of classical statistical mechanics 



II. KINETIC EQUATION 
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Let us consider the most general relativistic equation 
imposing the particle conservation during the evolution 
of a many body system. That equation is the first equa- 



2 



tion of the BBGKY hierarchy i.e. 



dl_ 

dp" 



d^p' (Ppi d^p[ 



(2.1) 

and describes, through the one-particle correlation func- 
tion or distribution function / = f{x,p), a relativistic 
particle system in presence of an external force field. The 
streaming term as well as the Lorentz invariant integra- 
tions in the collision integral, have the standard forms 
forms of the relativistic kinetic theory |46| . |47| . The two 
particle correlation function F{f,fi) [48!|, here denoted 
by f®fi, at the moment remains an unknown function, 
which will be determined in the following by employing 
the relativistic MCH. We recall that in classical kinet- 
ics, the two-particle correlation function, according to 
the MCH, is assumed to be, /(SD/i — f fi- 

Following standard lines of kinetic theory, we note that 
in stationary conditions, the collision integral in Eq. (|2.ip 
vanishes and then it follows 



(2.2) 



More in general it holds L{f®fi) = L{f'®f[), being L{x) 
an arbitrary function, and this relationship expresses a 
conservation law for the particle system. On the other 
hand a conservation law has form 



A(/) + A(/i)=.A(/') + A(/0 



(2.3) 



A(/) being the collision invariant of the system. There- 
fore we can pose L{f®Ji) = A(/) -I- A(/i). 

From the definition of the correlation function and tak- 
ing into account the indistinguishability of the particles, 
it follows that /<8)1 — 1 ® f = f and this implies that 
L(/) = A(/) and A(l) = 0. Therefore the function A(/) 
permit us to determine univocally the correlation func- 
tion through the equation 



A(/®/i)-A(/) + A(/i) 



(2.4) 



The composition law f <E) fi has the properties: 

(i) {f ® fi) ® f2 = f ® ifi "X" /2), i-e. it is associative, 

(ii) f<E)l = l<E)f = f, i.e. it admits as neutral element 
the unity, 

(iii) f ® fi = fi® f , i.e. it is commutative. 

The above three properties are the same of the ordi- 
nary product. In order to become ® a generalized prod- 
uct isomorphic to the ordinary product, it is necessary 
and sufficient to require the property 

(iv) f<E){l/f) = (1//)®/= 1, i.e. the inverse element of 
/ is 1//. Then, the real, positive, probability distribution 
functions, form an Abelian group. 

The function A(/), through Eq. (|2.4p . defines univo- 
cally the two-particle correlation function f ® fi which 
is a generalized product having the same algebraic prop- 
erties of the ordinary product. 

It is remarkable that the property (iv) imposes to the 
function A(/) the condition 



A(l//) = -A(/) 



(2.5) 



In relativistic kinetics, the collision invariant A(/), un- 
less an additive constant, is proportional to the micro- 
scopic relativistic invariant /. Then we can pose 



A(/) = -/3/ + /3Ai 



(2.6) 



/3 and being two arbitrary constants. In presence of 
an external electromagnetic field A", the more general 
microscopic relativistic invariant /, has a form propor- 
tional to 



/ = {p" + qA^/c) - mc^ , 



(2.7) 



t/y being the^ hydrodynamic four-vector velocity with 
VU^ = c2 |46|. 

Finally, after inversion of Eq. 
tionary distribution as follows 



, it obtains the sta- 



(2.8) 



It is remarkable to note that Eq. (|2.6p . after being 
written in the form — A(/) — /?/ + jSp, — 0, follows from 
the variational equation 
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S - 13 d'plf + f3fi d'pf 







(2.9) 



where the functional S, unless an arbitrary additive con- 
stant, is given by 



d'p aif) df , 



where the entropy density a{f) is given by 



A(/) = - 



df 



(2.10) 



(2.11) 



The variational equation (j2.9p . defining univocally 
the distribution (|2.8p . represents the Maximum En- 
tropy Principle (MEP). The constants (3 and are 
the Lagrange multipliers while the functional 5, defined 
through Eq. (j2.10p , is the system entropy. We stress that 
both the stationary distribution (|2.8p and the entropy 
(I2.10p depend on the form of the function A(/). This 
function is related to the two-particle correlation func- 
tion through ()2.4p . and its expression will be obtained in 
the next sections. 



III. THE H-THEOREM 

In the standard relativistic kinetics it is well known 
from the H-theorem that the production of entropy is 
never negative and in equilibrium conditions there is no 
entropy production. In the following we will demonstrate 
the H-theorem for the system governed by the kinetic 
equation (|2.ip . We define the four- vector entropy S'^ = 
{S°,S) as follows 



— P cr(/) , 



(3.1) 
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and note that 5*° = 5* coincides with the scalar entropy 
defined in Eq. (I2.10p . while S is the entropy flow. The 
identity (Pp/p^ — d^p 2 6'(p") 5{p^^p^ — m^c^), permits us 
to write the four- vector S'^ in the form 



5" 



d^p2 9ip")Sip'^p^-m^c^)p'' aif) . (3.2) 



In this expression d'^p is a scalar because the Jacobian of 
the Lorentz transformation is equal to unit. Then since 
transforms as a four-vector, we can conclude that S"^ 
transforms as a four- vector. 

In order to calculate the entropy production d^S'^ we 
start from the definition of S"^ given by (|3.ip . the evolu- 
tion equation (j2.ip and take into account the relationship 

aif) = [d<jif)/df]d,f = -A{f)d,f. We obtain 



d^p 

—FT InK(g) p" d^f 

f?^^^GI/«/!-/«/.lA(/) 



dp" 



(3.3) 



Since the Lorentz force has the properties p'^Fi, = 
and dF^ /dp^ = the last term in the above equation 
involving F'^ is equal to zero [4^. Given the particular 
symmetry of the non vanishing integral in Eq. p.3p we 
can write the entropy production as follows 



^^^^ G 



pQ p,0 



x[A(/) + A(/i 



Pi P'l 
-A(/') 



Hfi)] 



(3.4) 



Finally, we set this equation in the form 



" "47 p'O pf p'l> 

X [f'^f'i-f^fi] [A(/'®/0 - A(/®/i)] . (3.5) 

After imposing that A{h) is an increasing function, it 
results [hi — ^12] [A(/ii) — A(ft,2)] > V/ii, /12 and then we 
can conclude that 



d^S" > 



(3.6) 



This last relation is the local formulation of the relativis- 
tic H-theorem which represents the second law of the 
thermodynamics for the system governed by the evolu- 
tion equation (12. ip . 

In the next section we will determine the form of the 
function A(/) by employing the relativistic MCH. 



IV. RELATIVISTIC STATISTICAL THEORY 

In classical statistical mechanics, it is well known that 
the MCH, imposing for the two-particle correlation func- 
tion the classical form f ® fi = / /i , implies according 



to Eq. (|2.4[) . that A(/) = ln(/). As a consequence the 
expression of the entropy (|2.10p . unless an additive con- 
stant, simplifies as 



S ■ 



J d^pfA{f/e) 



(4.1) 



with 7=1 and e = e. It is worth stressing that the 
latter simple form of the entropy, defined as the mean 
value of — A(//e), is enforced exclusively by the MCH, 
and represents a very important feature of the ordinary 
statistical theory. 

In order to better point out the role of the MCH, in 
the construction of the ordinary statistical theory, we 
equal the expression (|2.10p . defining the entropy, with 
its expression as given in Eq. (|4.ip accounting the MCH, 
and obtain the following differential-functional equation 



-^/A(/) = i A(6/) 
df 1 



(4.2) 



The latter equation, imposed by the MCH, can be 
viewed as the starting point of an alternative approach 
to obtain the ordinary statistical theory. Indeed, as first 
step, we obtain the solution A(/) = ln(/), 7 = 1, e = e of 
Eq. (|4.2I) . obeying the condition (|2.5I) . As second step, 
through Eq. (|2.4p . we obtain the expression of the cor- 
relation function f®fi = f fi appearing in the evolution 
equation (|2.ip . As third step, by employing Eq. (|4.1I) . we 
obtain the Boltzmman entropy S = — J d'^p f ln(//e). 
As fourth and last step, through Eq. (j2.8p . we obtain the 
Maxwell distribution / = exp ( — /3/-I-/3/J,). 

It is important to emphasize that Eq. (|4.2p . already 
known in the literature, is enforced exclusively by the 
MCH. In the light of this, spontaneously emerges the 
question if Eq. (|4.2p with the condition (|2.5p admits 
other solutions, beside the aforementioned classical so- 
lution. In other words we pose the question if there 
exist other expressions of the A(/), apart the classical 
ones, preserving the important feature (|4.2p of the ordi- 
nary statistical mechanics, enforced by the MCH. Inter- 
estingly, Eqs. (14.21) and (12. 5p admit another, more gen- 
eral solution, which is unique and contains as limiting, 
particular case, the classical solution. The prediction of 
this new solution represents a very strong argument in 
favor of the present approach. 

The general solution of Eq. (|4.2p is defined up to an 
arbitrary multiplicative constant which can be fixed by 
using the normalization condition A'(l) = 1 yielding 



7 = A(e) 



(4.3) 



Consequently the general solution, depends only on one 
free parameter, in the following indicated with k, and 
defines a generalized logarithm i.e. 



(4.4) 



After solving Ea. (|4.2l) with the condition (|2.5p . one ob- 
tains [2] the following simple expressions for the general- 
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ized logarithm and its inverse function i.e. the general- 
ized exponential A~^{x) = exp^(a:): 



which in the k — ^ classical limit, reduces to ffi. 



2k 



.(x) = + k'^X^ + KX^ , 



exp, 



(4.5) 
(4.6) 



with < 1. The constants 7 and e, depend on the 
parameter k and are given by 



1 
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e = exp^(7). 



(4.7) 



The above generalized logarithm and exponential are 
one-parameter continuous deformations of the ordinary 
functions which recovers in the k — limit. 

The entropy (|4.1I) associated to the new solution be- 
comes 



(4.8) 



S = -jj d^p /ln,(//e) , 
and can be written explicitly as 

while the stationary distribution (12. 8p assumes the form 

/ = exp,(-/3/-f . (4.10) 

The latter entropy and distribution, in the k — > limit 
reduce to the classical Boltzmann entropy and Maxwell- 
Boltzmann distribution respectively. 

The distribution (|4.10l) in the global rest frame where 
Uv = (c, 0,0, 0) and in absence of external forces i.e. 
A'^ = 0, simplifies as 



/ = exp, i-PE + Pl,) 



(4.11) 



E being the relativistic kinetic energy. This distribution 
at low energies {E 0) reduces to the classical Maxwell- 
Boltzmann distribution i.e. / « exp {—/3E — /?/i), while 
at relativistic energies (E — ^ -l-c») presents power-law 
tails / oc E~^^'^, in accordance with the experimental 
evidence. 

After posing A(/) = lnK(/), in Eq. (|2.4|) . the relativis- 
tic two-particle correlation function assumes the form 



/ » /i = exp^ ( In^/ + In^/i ) 



(4.12) 



It is easy to verify that the distribution (I4.10p is the 
stationary solution of the evolution equation (|2.1|) . First 
we observe that in stationary conditions, both the steam- 
ing term (i.e. the left hand side) as well as the collision 
integral (i.e. the right hand side) of Eq. ()2.1|) . vanish. 

For the collision integral to vanish, it is necessary and 
sufRcient to pose f ® f'l — f® fi — 0. The latter equation, 
after taking into account Eqs. (|4.1ip and (|4.12l) . produces 
the relationships E' + E[ = E + Ei and /.i' + fJ.'i = /i -f /ii , 
expressing the particle energy conservation and the par- 
ticle number conservation respectively, during the colli- 
sions. 

On the other hand, from relativistic kinetic theory |46| . 
it is well known that the streaming term, in the evolu- 
tion equation, vanishes when the distribution function 
depends on the microscopic relativistic invariant, given 
by (j2.7|) . which in the case of the distribution (|4.11|) re- 
duced to I = E. 

We recall that the first experimental validation of the 
distribution (|4.1ip . concerns cosmic rays and has been 
considered in ref. Recently a computer validation 

of the same distribution, has been considered in refs. 
[20I [2l| , where the relaxation in relativistic plasmas un- 
der wave-particle interaction, has been simulated numer- 
ically. 

As conclusions we emphasize that: 

(i) The results obtained in sections II and III hold for 
a very large class of relativistic statistical theories, in- 
cluded the Juttner theory. The arbitrary function A(/), 
after being fixed, defines univocally the statistical the- 
ory. Indeed, after fixing the form of the function A(/), 
the two-particle correlation function defined through Eq. 
(|2.4p . the evolution equation (|2.ip . the stationary dis- 
tribution function (j2.8p . and the four- vector entropy de- 
fined through Eqs. (|2.1ip and (13. ip . are univocally fixed. 
Moreover during the evolution of the system the second 
law of the thermodynamics holds at any time, indepen- 
dently on the particular form of the function A(/). 

(ii) In fourth section, the general results obtained in 
sections II and III, are employed in order to construct 
the kinetic foundations of the statistical theory predicting 
the distribution (|4.1ip . 
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